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Abstract 

In this paper, taking the view that a linear parallel interference canceller (LPIC) can 
be seen as a linear matrix filter, we propose new linear matrix filters that can result 
in improved bit error performance compared to other LPICs in the literature. The 
motivation for the proposed filters arises from the possibility of avoiding the generation 
of certain interference and noise terms in a given stage that would have been present 
in a conventional LPIC (CLPIC). In the proposed filters, we achieve such avoidance 
of the generation of interference and noise terms in a given stage by simply making 
the diagonal elements of a certain matrix in that stage equal to zero. Hence, the 
proposed filters do not require additional complexity compared to the CLPIC, and 
they can allow achieving a certain error performance using fewer LPIC stages. We 
also extend the proposed matrix filter solutions to a multicarrier DS-CDMA system, 
where we consider two types of receivers. In one receiver (referred to as Type-I receiver), 
LPIC is performed on each subcarrier first, followed by multicarrier combining (MCC). 
In the other receiver (called Type-II receiver), MCC is performed first, followed by 
LPIC. We show that in both Type-I and Type-II receivers, the proposed matrix filters 
outperform other matrix filters. Also, Type-II receiver performs better than Type-I 
receiver because of enhanced accuracy of the interference estimates achieved due to 
frequency diversity offered by MCC. 
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detector, multicarrier DS-CDMA. 
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1 Introduction 

Linear parallel interference cancellers (LPIC) have the advantages of implementation simplic- 
ity, analytical tractability, and good performance p]-[12]. The conventional way to realize 
LPIC schemes is to use unsealed values of the previous stage soft outputs of different users 
for multiple access interference (MAI) estimation. In [3], Guo et al described and analyzed 
LPIC schemes for CDMA using a matrix-algebraic approach. They pointed out that an 
LPIC can be viewed as a linear matrix filter applied directly to the chip matched filter (MF) 
output vector. While the matrix filter corresponding to the conventional LPIC (CLPIC) con- 
verges to the decorrelating (DC) detector, they also proposed a modified matrix filter which 
converges to a minimum mean square error (MMSE) detector. This was done by exploiting 
the equivalence of the LPIC to a steepest descent optimization method for minimizing the 
mean square error. For this optimization, they obtained optimum step sizes for different 
stages that remove the excess mean square error in K stages (where K is the number of 
users), leaving only the minimum MSE in stages greater than K. The condition for this 
convergence has been shown to be that the maximum eigenvalue of the correlation matrix 
must be less than two. 

Our contribution in this paper is that we propose new linear matrix filters that can perform 
better than the matrix filters studied in [3]. The motivation for the proposed filters arises 
from the possibility of avoiding the generation of certain interference and noise terms in a 
given stage that would have been present in the CLPIC. In the proposed filters, we achieve 
such avoidance of the generation of interference and noise terms in a given stage by simply 
making the diagonal elements of a certain matrix in that stage equal to zero. Hence the 
proposed filters do not require additional complexity compared to the CLPIC. We show that 
the proposed matrix filters can achieve better performance compared to other matrix filters 
in the literature. This, in turn, can allow achieving a certain error performance using fewer 
LPIC stages. We also propose filters that use different step sizes for different stages (but the 
same step size for all users at a given stage). In addition, we propose filters that use different 
weights for different users in different stages, where we also obtain closed-form expressions 
for the optimum weights that maximize the output average SINR in a given stage. 
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We further extend the proposed matrix filter solutions to a multicarrier DS-CDMA system, 
where multicarrier combining (MCC) needs to be carried out in addition to the LPIC op- 
eration. Depending on which operation (i.e., MCC or LPIC) gets done first at the receiver, 
the resulting performances and complexities differ. We consider two types of receivers. In 
what we call the Type-I receiver, LPIC is performed on each subcarrier first, followed by 
MCC [12] . In the Type-II receiver, MCC is performed first, followed by LPIC. We show that 
in both Type-I and Type-II receivers, the proposed matrix filters outperform other matrix 
filters. Also, Type-II receiver outperforms Type-I receiver because of enhanced accuracy of 
the MAI estimates achieved due to frequency diversity offered by MCC. 

The rest of the paper is organized as follows. In Sec. [2, we present the system model. In 
Sec. El we present the proposed matrix filters for single carrier DS-CDMA, along with their 
bit error performance results. Section H] presents the proposed filters and their performance 
in multicarrier DS-CDMA. Conclusions are presented in Sec. [5j 



2 System Model 

We consider a .fT-user synchronous multicarrier DS-CDMA system with M subcarriers. Let 
bk G {+1, —1} denote the binary data symbol of the kth user, which is sent in parallel on 
M subcarriers [13]. Let P denote the number of chips-per-bit in the signature waveforms. 
It is assumed that the channel is frequency non-selective on each subcarrier and the fading 
is slow (assumed constant over one bit interval) and independent from one subcarrier to the 
other. We assume that the subcarriers are separated sufficiently apart so that inter-carrier 
interference is negligible. 

Let yW'W = [Vi 1 ^*'* y-p'® ■ ■ ■ y^'^ denote the i^-length received signal vectoJ3 at the 



MF output on the ith subcarrier; i.e., y^ hW is the MF output (i.e., 1st stage output) of the 
kth user on the ith subcarrier, given by 



1 Vectors are denoted by boldface lowercase letters, and matrices are denoted by boldface uppercase letters. 
[.} T and [.] H denote transpose and conjugate transpose operations, respectively. Re{a} and Im{a} denote 
the real and imaginary parts of a. 

2 We take the MF output (i.e., the despread output) as the 1st stage output in the multistage LPIC 
receivers. So, the ( m )'W in the superscript of y denotes the stage index m and subcarrier index i. 
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The MF output vector yW>W can be written in the form 



R»H»b + n« 



where represents the KxK channel matrix, given by 
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and the channel coefficients ft, 
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where = p?/ is the normalized cross-correlation coefficient between the signature wave- 
forms of the Ith and jth users on the ith subcarrier. The iT-length data vector b is 

T 



given by b 



A^A 



where Ak denotes the transmit amplitude of 



the kth user. The l"f -length noise vector nW on the ith subcarrier is given by n 



,0 



/■?. 



A 



, where ni 1 "* denotes the additive noise component of the kth user 



on the ith subcarrier, which is assumed to be complex Gaussian with zero mean with 



a 2 when j = k and E[n 



in 



a2 Pkj wnen j k. 



3 Proposed Linear Matrix Filters 

In this section, we propose improved LPICs for single carrier DS-CDMA (i.e., M — 1 in 
the system model presented in Sec. EJ. So, for notational simplicity, we drop the subcarrier 
index (i) in this section. In Sec. HJ we will consider M > 1. We assume that all the channel 
coefficients are perfectly known at the receiver. Dropping the subcarrier index in (JTJ), the 
MF output of the desired user k, y[ , can be written as 
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Vk - x k + x jPkj+ n>k , where xi = Afyhi, l = l,2,---,K. (5) 

desired signal j t noise 

MAI 

3.1 Conventional Matrix Filter, G (m) 

In conventional LPIC (CLPIC), an estimate of the MAI for the desired user in the current 
stage is obtained using all the other users' soft outputs from the previous stage for cancella- 
tion in the current stage. The mth stage output of the desired user k, y^, in CLPIC is [9] 

vt ] = vP- E p ]k y ( r 1] - (6) 

j=i,j¥=k 

s * ' 

MAI estimate 

The fcth user's bit decision after MAI cancellation in the mth stage, b^\ is obtained as 

6i m) = sgn(Re(^i m) )). (7) 
The CLPIC output in ([6]) can be written in matrix algebraic form as [3] 

m 

y « = £(i-Rr x y (1) ; (8) 

3=1 

~ - v ' 

G( m ) 

The G (m) filter in ([8]) can be viewed as an equivalent one-shot linear matrix filter for the 
mth stage of the CLPIC. 

3.2 Proposed Matrix Filter, G p (m) 

In this subsection, we propose a new linear matrix filter, which we denote as Gp^, that 
can perform better than the matrix filter G*" 1 -' in (jSj). The motivation for the new matrix 
filter can be explained as follows. 

What does the matrix filter do: It is noted that the behavior of the G^™' filter in (jHJ) 
(i.e., CLPIC) at a given stage m > 2 is characterized by a) interference removal, b) generation 
of new interference terms, c) desired signal loss/gain, d) desired signal recovery/removal, and 
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I + (I - R) + (I - R) 2 + • ■ ■ + (I - R 
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e) noise enhancement. For example, the cancellation operation in the 2nd stage (i.e., m — 2) 
results in i) interference removal, ii) generation of new interference terms, Hi) desired signal 
loss, and iv) noise enhancement. This can be seen by observing the 2nd stage output 
expression for the desired user k, which can be written, using and ([5]), as 



(2) 

Vk 
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desired signal loss new interference terms additional noise terms 

Comparing the expression at the MF output, in (jSJ) and the expression for the 2nd stage 
output, \ in ([9}, it can be seen that the cancellation operation in the 2nd stage results in 
the following at the 2nd stage output. 

• The interference terms, J2j^k PjkXj, that were present in the MF output in ((SJ) are re- 
moved. In the process, i) new interference terms proportional to p 2 , i.e., J2j^ k Pjk Pji x i 
in (jUj), get generated, ii) a fraction J2j^kP]k °f the desired signal component gets lost, 
and Hi) additional noise terms proportional to p, i.e., J2j^kPjk n j i n ©, get introduced. 

In Appendix A, we present the expression for the 3rd stage output in an expanded form. 
From (|4ip in Appendix A, we can make the following observations which result from the 
cancellation operation in the 3rd stage. 

• The desired signal loss that occurred in the 2nd stage is recovered (see the two [T] terms 
cancelling each other in (|4ip ). In the process, new interference terms proportional to 
p 3 (see the | bj | term in (1411) ) as well as additional noise terms proportional to p 2 (see 
the | b n | term in (I4T1) ) get generated. 

• Interference terms generated in the 2nd stage are removed (see the two [c] terms can- 
celling each other in (14T1) ). In the process, i) further desired signal loss/gai propor- 



3 Depending on p's being positive or negative, the term [~d] in (|¥T|l can be positive or negative, because 
of which there can be a desired signal gain or loss. 



tional to p 3 occurs (see the [d] term in (|4ip ). and ii) new interference terms proportional 
to p 3 (see the | ej | term in (|4T1) ) as well as additional noise terms proportional to p 2 
(see the | e n | term in fid"]) ) get generated. 

Similar observations can be made on the expanded form of the equations for the subsequent 
stages of the CLPICfj. For m — > oo, the CLPIC is known to converge to the decorrelating 
detector, provided the eigenvalues of the R matrix are less than two |3J. That is, when 
m — > oo, in the expression for G^ output in ( l4"5j) . the desired signal loss/gain and the 
interference terms go to zero and the noise term gets enhanced. 

What is proposed to be achieved using the G p ^: As explained above, in the 
filter, new interference and noise terms get generated in the process of interference removal 
and recovery/removal of desired signal loss/gain. We seek to avoid the generation of some of 
these new interference and noise terms. For example, as will be shown next, the generation 
of the | bj | and | b n | terms at the 3rd stage output in (14"T1) can be avoided by simply making 
the diagonal elements of a certain matrix in the cancellation operation in the 3rd stage equal 
to zero. This, as we will see later, can result in improved performance compared to the G^ 
filter. 

Proposed matrix filter, G p ^ m ^: We propose to avoid the generation of new interference 
and noise terms in T 3 in (jSJ), caused in the process of recovery /removal of desired signal 
loss/gain in the previous stage. Since there is no desired signal loss/gain in the 1st stage, the 
2nd stage of the proposed filter is the same as that of the G^ filter, i.e., Gp*- 2 - 1 = G^ 2 -*. For 
stages greater than two, i.e., for m > 3, the mth stage output of the proposed filter G p ^ m \ 
denoted by , can be written as 

vfi? = vti 1] + (-!) m+1 EE E ■•■ 

kiy^k k2^k,k\ kgy^k,k2 

K K 

E E ■ ■ ■ Pk 3 k 2 Pk 2 k 1 yk 1 ■ (10) 

We note that the above expression is obtained by i) dropping T3 from (jSj), and ii) modifying 
T4 in ()44l) such that all the summations in it exclude the desired user index k. The above 



4 The general expression for the mth stage output in expanded form, for any m > 3, and the corresponding 
observations are given in Appendix B. 



two modifications ensure that the proposed filter removes the previous stage interference 
while avoiding the recovery/removal of the desired signal loss/gain^. Also, because of these 
modifications, the interference and noise terms in a given stage of the proposed filter will be 
a subset of the interference and noise terms in the same stage of the filter. Equation 
([TO]) can be written in the following form 

K / K K K 

VkZ> = Uk ~ E I P fcfc i ~ E Pkk 2 Pk 2 k 1 + E E Pkk 3 Pk 3 k 2 Pk2ki — ••• 
k\^k \ k 2 ^=k,k\ k 2 ^k,ki k 3 ^k,k 2 

K K K K \ 

+ (-1)™ E E "' E E Pkk m - X Pk m -ik m - 2 ■ ■ ■ Pk 3 k 2 Pk 2 kx) Vkh 

k 2 ^k,k\ k^k,k 2 k m - 2 ^k, 

which, in turn, can be expressed in matrix form as 

(m-l \ 
EBij y {1) , (12) 



G p ( ra ' 



where 

B 



B^CE-R)! , (13) 



[M] denotes the matrix M with its diagonal elements made equal to zero, and B = I. 
Note that, since ( fl2i) is structurally the same as (JED and the [.]° operation in ( 1T3|) does not 
add to complexity, the proposed G p ' m ' filter has the same complexity as the G^ filter. 

The G^" 1 ) filter is known to converge to the decorrelating detector for m — > oo, provided 
the maximum eigenvalue of the R matrix is less than two [3]. That is, G^ = R -1 , which 
results in the output vector 

(y(°°>) G = R-V (1) = x + R^n. (14) 

As with G^ m \ all the interference terms in G p ^ also go to zero for m — > oo. This can be 
seen as follows. From (fl2|) and ( [TBI . Gp^ 00 - 1 can be written in the form 

G p (00) = I +[(I-R)-D 1 ] + {[(I-R)-D 1 ](I-R)-D 2 } + ... J (15) 

B ° Bi B 2 



5 Although possible signal loss recovery is avoided in the process, the net effect can still be beneficial (we 
will see this in Sec. 
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where D n is a diagonal matrix with the diagonal elements the same as those in the matrix 
B n _i(I — R). Equation f|T5|) can be written as 

G p (°°) = (I + (I - R) + (I - R) 2 + • • •) - D x (I + (I - R) + (I - R) 2 + • • •) 

R 1 



(16) 



- D 2 (I+(I-R) + (I-R) 2 + ---) - 
(I-Da-D 2 - • •) R 1 . 



Hence, the output vector for m — > oo is given by 

(yM) = FR-V (1) = Fx + FR ^n. 

V / G p 

The diagonal matrix F defined in (fTBTl can be written as 



diag\Jx,f 2 i---jK 



where is given by 



(17) 



18) 



K K K K K K 

fk = 1 - X! Pkk x Pk x k + ^2 X! Pkk 2 Pk 2 k l Pk 1 k - ^2 X! ^2 Pkk 3 Pk3k 2 Pk2k 1 Pkik ' ' ' (19) 
ki^k k\^kk 2 -^k,k\ kij^k k 2 j^k,ki kg^k,k 2 

For the case of equi-correlated users, fk in f[T§j) can be shown to converge to 

1 — ((K — l)p 2 /(l + (K — 2)p), and there are no interference terms in (fTTj) . Also, note that 

the outputs of the G filter in (11411 and the G p filter in (TTTj) have the same SNR for m — ► oo. 



3.3 Why G p (m) can perform better than - An Illustration 

To analytically see why G p (m ' ) can perform better than G^ m ^ for m > 3, consider the case 
of K equi-correlated users with correlation coefficient p, and no noise. Let us consider the 
average signal-to-interference ratio (SIR) at the 3rd stage output for G p ^ m ' and G*"^. The 
3rd stage output of G^ m \ in the absence of noise, can be written as 



1 + 



D 



Xk + 



Bi 



+ 



(20) 



where the terms 
can be written as 



D 



Bi 



Hi 



are defined in (14TI) . Likewise, the 3rd stage output of G p < - m ' ) 



(3) 

Vh 



+ 



D 



X k + 



(21) 



Note that the interference term I Bj I generated in G*" 1 -' is not generated in G p ^. Also, the 
desired signal term [X] is recovered in G ( - m - ) whereas it is not recovered in G p ^ m \ Now, from 
(120]) . the average SIR at the 3rd stage output of G^ m \ for the case of equi-correlated and 



equal-amplitude users, can be obtained as 



(3) l + (A-l)(A-2)p 



2 



SIRL= L T2- (22) 



G 



(A-l)((A-l)p 3 + (A-2) 2 p 3 , 
Likewise, from (12ip . the average SIR at the 3rd stage output of G^™^ can be obtained as 



(AT- l)((A-2) 2 p 3 
From (1221) and (|23|) . it can be seen that 



( SIR ) Gp = 1 | (A-l)(l-(A-2)p)(l + (A-2)p-(A-l)p 2 ) 
J (STR)^ + - 2) 2 (l + (K-1)(K- 2)/ 



(24) 



For p > 0, the maximum eigenvalue for the R matrix is 1 + (A — l)p, so that the condition 
for convergence is 1 + (A — l)p < 2, i.e., (A — l)p < 1. Now, in ( 1241) . the 2nd term is 
positive when (A — l)p < 1, which results in /3 > 1. This implies that G p ( - ,?1 ' ) results in a 
higher average output SIR than G^ m \ As will be seen in Sec. I3.6[ simulation results show 
that Gp^ can perform better than G^ m ) in the case of non-equicorrelated users as well as 
in the presence of noise. 



3.4 A Modified MMSE Converging Filter, G 



(m) 



As pointed out in Sec. [fl Guo et al, in [3], have proposed modifications to the G^ filter 
so that the resulting modified matrix filter converges to the MMSE detector instead of the 
decorrelating detector, by exploiting the equivalence of the LPIC to the steepest descent 
method (SDM) of optimization for minimizing the MSE. They also derived optimum step 
sizes for various stages, which ensured convergence to the MMSE detector in A stages, where 
A is the number of users. We refer to this MMSE converging matrix filter proposed by Guo 
et al in [3] as the G^ 1 filter, which is given by [3] 
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(1) 



(25) 



1=1 



i=i 



where \i{ is the step size at stage i, the optimum values of which were obtained to be 

1 



,. . (26) 
Xi + a 2 

where Aj, % — 1, 2, • • • , K are eigenvalues of matrix R. We note that a similar SDM view can 
be taken to modify our proposed matrix filter G p ^ so that it can converge to the MMSE 
detector. We refer to such a modified version of our proposed filter as G p ^ m ^ filter, where 
we avoid the generation of new interference and noise terms as in G p ^ m \ while using the step 
sizes obtained for G^ m ^ in [3J. Accordingly, we propose the modified version of the MMSE 
converging filter as 



jpfi 



m—l 



l^m—i"i I y 



(1) 



(27) 



G(m) 
PM 



where J« is given by 



Ji_l I - /! A '-i+l(R + O I 



and J = I. 



(2f 



3.5 A Weighted Matrix Filter, G p ^ m) 

In the G^ m ' and G p ^ m ^ filters above, different step sizes are used in different stages (but the 
same step size for all users in a stage). Improved performance can be achieved if different 
scaling factors (weights) are used for different users in different stages. Accordingly, we 
propose a weighted version of our proposed filter G p ^ m \ We refer to this weighted version 
as G pw ^ m \ and is derived as follows. 

In a weighted LPIC (WLPIC), the MAI estimate in a given stage is scaled by a weight before 
cancellation (unit weight corresponds to CLPIC and zero weight corresponds to MF). For 
example, the mth stage output of the desired user k, yj^, in a WLPIC is given by 

K 



(m) 



(1) (m) 

Vk ~ w k 



PikVj, 



(m-l) 

w ) 



(29) 
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where w k m ^ is the weight with which the MAI estimate for the kth user in the mth stage is 
scaled. For m > 2, the weighted cancellation operation in ( |29i) can be written in the form 

y£"> = (l + W (m) (I - R) + W (m) (I - R)W (m_1) (I — R) H 

+ W (m) (I - R) W (m - 1} (I - R) • • • W (2) (I - R)^ y (1) , (30) 

where W( m ) is the weight matrix at the mth stage, given by W< m ) = diag ( 



w 2 ), 

and = 0. Now, as in G^ m \ in order to avoid the generation of new interference and 

noise terms, we modify (13"U1) as follows: 

yS = I E b, ] y (1) , (3i) 



where 

a 




;B n _!)(W( m - n+1 ))(I - R)l° , and B = I. (32) 



Note that G pw (m) becomes G p (m) when W( m ) = I, Vm > 1. 

Optimum Weight Matrix, W^: The mth stage output of the kth user when G pw < - m * ) 
filter is used can be written as 

K 

( (m)\ (1) (m)V^ (™) (!) 

[Vk L = Vk w k 2^Qk,iVi 
Gpw i^k 

( K \ K K ( ( K \\ 

I i (™) ( m ) \ , ( m ) (™) , / ( m ) / ( m ) , ( m ) I I /oo\ 

= a>fcll-ti;£ 2^ q k,i Pki)+ n k -w\ '2^1k/ n ^ + l^\Pki-wl '\ql/+ }^ 8*,fc>*i* I Jar*, (33) 

\ i^k J i^Lk i^k \ V kijii,k / / 

v ' S v ' V v ' 

desired signal noise interference 

where 

q(™) _ ^p k - ^ -,.,( m_1 )^.. . _i_ ,„( m - 2 ) .„( m_1 ). 
if 



Since the interference and noise terms in on the RHS of f)33p are the sum of linear combi- 
nations of complex Gaussian r.v's (since the fade coefficients hk are assumed to be complex 
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Gaussian), the average SINR for the kth user at the mth stage output can be obtained, in 
closed-form, as 



awl 



where 



(r„\ ALII LbU I. 

SINR™ = -^- 2 (35) 



K K K K 

i^A; i^fc i^fc k\^i,k 

K K K K 

i^fc k\^i,k ij^k jj^k 



By differentiating the average SINR expression in fl35l) w.r.t. and equating to zero, the 
optimum weights w^ pt can be obtained, in closed-form, as 



™W _ rf ~ a6 ( -> )()) 

"W - c _ ad + a 2 (e _ a 2Y W 



In Fig. [T], we plot the average SINR at the mth stage output of the proposed G pw (m ' ) filter, 
as a function of weight, w^, at an average SNR of 20 dB. The average SNR of user k is 
defined as A\ja 2 . The number of users considered is K = 20, the processing gain is P = 64, 
and there is no near-far effect (i.e., A\ — A2 — ■ ■ ■ = Art). From Fig. [TJ it can be seen 
that for a given stage index m, the maximum output average SINR occurs at an optimum 
weight; the closed-form expression for this optimum weight is given by fl36|) . The maximum 
average SINR increases as m is increased. Also, we see diminishing improvement in SINR 
with increasing m, as expected. Another key observation in Fig. [T]is that, while non-unity 
weights are optimum for small values of m, the optimum weights approach unity for large m. 
G pw being structurally similar to G p except for the weights, like G p and G filters, G pw is 
also expected to converge to R 1 for m — > 00, and this explains why w^ pt — > 1 for m — ► 00. 

3.6 Results and Discussion 

In this subsection, we present a comparison of the bit error rate (BER) performance of 
different matrix filters. The various matrix filters considered include: i) the conventional 
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filter, G( m \ given by ©, ii) the proposed filter, G p (m) , given by (JT2D, Hi) the MMSE 
converging filter in [3], G^ m \ given by fl25|) . if) the modified MMSE converging filter, 
G p/ / , given by ( 1271) . and v) the proposed weighted filter, G pvf ^ m \ given by ( 13~TT) . 



In Fig. [21 we plot the BER performance of the conventional filter, G^ m \ and the proposed 
filter, G p ( m \ as a function of the stage index, m, for M — 1, K — 20, P = 64, and average 
SNR = 15 dB, for both no near-far (i.e., Ax — Ai — ■ ■ ■ — Ak) as well as near-far conditions. 
In all the simulations, user 1 is taken to be the desired user. Random binary sequences are 
used as spreading sequences. For the near-far condition, odd-indexed users (users 3, 5, 7, • • •) 
transmit with the same amplitude as the desired user 1, whereas the even- indexed users (users 
2, 4, 6, • ■ ■) transmit at 10 times larger amplitude than the desired user. The performance of 
the MF detector and the DC detector are also plotted for comparison. From Fig. [2, it can 
be seen that the conventional G^ m ) filter approaches the DC detector performance rather 
slowly for increasing m. Observe that the performance of the proposed G p ^ filter and the 
conventional G^ filter are the same for m = 2 because of no desired signal loss recovery at 
the 2nd stage of both G^ and G p ^ m \ However, for m > 3, the G p ' m ' filter performs better 
than the G (m) filter. This is because the G p ^ filter, as intended, avoids the generation 
of new interference and noise terms (e.g., I b t I and I b n I terms for m = 3) compared to the 



q("») filter. The G p ^ filter is found to offer greater advantage in near-far conditions, since 
strong other-user interference terms in I Bj I are avoided in the G p ^ m \ 



Next, in Fig. [3], we present a comparison of the performance of the MMSE converging G At (m ' ) 
filter in [3], and the modified MMSE converging, G p ^ m \ for the same system conditions in 
Fig. [2l The performance of the MF and MMSE detectors are also plotted for comparison. 
Here again, the G^" 1 - 1 and G p/1 ' m ' filters perform the same for m = 2. Also, both G^" 1 - 1 and 
G p/1 ' m ' are seen to approach the MMSE performance as m is increased. For m > 3, G PAt (m ' ) 
performs better than G^ m ^ because of the avoidance of new interference and noise terms. 
In generating the plot for G PfJ- ^ m \ we have used the step sizes in f )26p . which are actually 
optimum for G^/"^. Even with these step sizes (which can be suboptimum for Gp^™-*), the 
proposed G p ^ m ^ filter approaches the MMSE performance faster than the G^" 1 -* filter. 

Finally, in Fig. HI we illustrate the performance of all the matrix filters considered in this 



14 



paper, including the proposed weighted filter, G pv / m % under the no near-far condition. The 
performance of the MF, DC and MMSE detectors are also plotted. It can be observed that 
among all the filters considered, the proposed weighted filter G pvr ^ performs the best for 
small values of m (m < 6, for example). In other words, G pw ^ performs best in terms 
of convergence, i.e., fewer stages are sufficient to yield close to DC detector performance. 
This may be expected, because in the G^™** and G p/ / m -' filters the optimum step sizes are 
obtained only on a per-stage basis, whereas in the Gpw 1 -"^ filter the optimum weights are 
obtained on a per-stage as well as a per-user basis. The computation of the optimum weights, 
w kopti f° r the Gpw 1 -"^ filter, using the closed- form expression in (1361) . adds to the receiver 
complexity. However, since these optimum weights are computed by using the average SINR 
expression, the weights computation can be carried out off-line once (or whenever users exit 
from or enter into the system, which changes the correlation matrix), and this need not 
add to the per-bit complexity of the canceller. In Fig. HI we also show the performance of 
the conventional weighted LPIC given in [9], denoted by G w (m) , for up to m = 4. As can 
be seen, because of the SIR maximization using optimum weights, the performance of the 
G w ' m ' filter is almost the same as that of the proposed G pw ' m ' filter. We further note that 
the optimum weights expressions for the G w ^ m ^ filter need to be derived separately on a 
stage by stage basis - the optimum weights expressions for up to m = 4 are given in [9],[10j, 
and the optimum weights derivation becomes increasingly cumbersome for increasing m. On 
the other hand, the feature of making the diagonal elements zero in the G pw (m) filter allows 
optimum weights expressions to be obtained for any m (given by Eqn. (1361) ). In terms of 
convergence as well as complexity, the proposed filter G PA / m ^ is also quite attractive. 

4 Proposed Filters in Multicarrier DS-CDMA 

In this section, we extend the proposed matrix filter solutions in the previous section to 
multicarrier DS-CDMA (i.e., for M > 2). Here, the multicarrier combining (MCC) operation 
has to be performed in addition to LPIC. Depending on which operation (i.e., MCC or LPIC) 
gets done first at the receiver, the resulting performances and complexities differ. 
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Type-I Receiver: We first consider a receiver where we perform LPIC first on each subcar- 
rier, followed by MCC as shown in Fig. [5] [12]. We refer to this receiver as Type-I receiver. 
We note that with this Type-I receiver using the conventional filter G^ m \ tractable BER 
analysis becomes feasible (reported in [12J). Also, this receiver architecture can be viewed as 
a direct adoption of the filters proposed for single carrier DS-CDMA in Sec. |3j on individual 
subcarriers in the MC DS-CDMA system. Hence, all the matrix filters in Sec. |3j namely, 
G p ( m \ G^\ G PA / m ^ and G pvv / m \ can be directly employed on the individual subcarriers. 

Type-II Receiver: Since MCC operation can provide frequency diversity, performing MCC 
before LPIC can enhance the accuracy of the estimates of the MAI and hence improve 
performance. Accordingly, we propose a Type-II receiver, where MCC is performed first, 
followed by LPIC, as shown in Fig. El The output of the MC combiner in vector form, 
denoted by y^ 1 -*, can be written as 

y c ^ = R c b + z, (37) 

where y^ = {y^ y} l) ■■■ yp] T , yf ) = y { ^\ (H»f R.W H « 

and z = \Y-U {hx®)*nf (h2®)*r$ ■■■ E^i (h K ^X n^Y . The conventional filter 



n K 



at the MCC output (referred to as G c ^ filter) has its mth stage output vector given by 

/ m— 1 \ 

= i+E(i-W y c(1) , (3« 



i=i 



where R eff = R C H D _1 , and H D = E 4 =i (n {i) ) H H {i) ■ Similar to the proposed G p (m) filter 
for single carrier DS-CDMA in the previous section (where the idea of zeroing the diagonal 
elements of a certain matrix is adopted), a proposed matrix filter for MC DS-CDMA Type-II 
receiver, denoted by G p c( - m \ can be obtained as 

y cH = I^B^y^), (39) 




where BJ; 



B^_! (I - R e ff) 



and B c = I. 



Results and Discussions: In Fig. [TJ we present the simulated BER performance com- 
parison between the Type-I and Type-II receivers for MC DS-CDMA using different matrix 
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filters. Specifically, we compare the performance of i) Type-I receiver with the G(' m ) filter, 
ii) Type-I receiver with the G p ( - m ' ) filter, Hi) Type-II receiver with the G^" 1 -* filter, and 
iv) Type-II receiver with the G p c< -"^ filter. We also compare the performance of the above 
detectors with DC and MMSE detectors. Random binary sequences of length P are used 
as the spreading sequences on each subcarrier, and the average SNR of user k is defined as 
— 2 K . In Fig. [71 we plot the BER as a function of stage index, m, for M — 4, K — 20, 
P = 64, and average SNR = 14 dB. The following observations can be made from Fig. [7J 

• Comparing the performance of Type-I and Type-II receivers for a given filter, we 
observe that Type-II receivers perform significantly better than Type-I receivers. For 
example, comparing the performance of Type-I receiver with the G^ filter and Type- 
II receiver with the G c< - m ) filter, we see that Type-II receiver with the G^" 1 -* filter 
performs significantly better (e.g., for m = 4, Type-I with the G^ filter gives a BER 
of 8 x 1(T 2 , whereas Type-II with the G c(m) filter results in a BER of 2 x 1(T 4 ). The 
superiority of Type-II receivers is consistent across all filters considered, i.e., G, G p , 
DC, MMSE. This superiority of Type-II receivers is expected, since the MAI estimates 
can be more accurate in Type-II, because of multicarrier combining before IC. 

• Like in SC DS-CDMA, in MC DS-CDMA, the proposed G p filter performs better than 
the G filter. This is observed to be true in both Type-I as well as Type-II receivers. 
For example, Type-I with the G p *- m ^ filter achieves a BER of 9 x 10~ 3 in just 5 stages, 
whereas the same BER is achieved by Type-I with the G ( - m - ) filter only after 15 stages. 
Likewise, in Type-II receivers, the G p c ^ filter performs better than the G c ^ filter. 

5 Conclusions 

We proposed improved LPICs for CDMA by viewing an LPIC as a linear matrix filter. 
Specifically, we proposed new linear matrix filters which achieved better performance than 
other linear matrix filters in the literature. This was made possible by avoiding the gen- 
eration of certain new interference and noise terms by making the diagonal elements of a 
certain matrix equal to zero in each stage, without adding complexity. We also extended the 
proposed matrix filter solutions to multicarrier DS-CDMA, where we considered two types 
of receivers; in both types of receivers the proposed filters were shown to outperform other 
filters in the literature. 
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Appendix A: Expression for the 3rd stage G (m) filter output 

In this appendix, we write the expression for the 3rd stage output of the G^ m ' filter (i.e., 
CLPIC) in an expanded form. From (JS|), can be written as 



r(3) 



I + (I - R 



,(2) 



(I-R) 2 y (1) 



y w = y v 

K K K K K 



= Vk +Y Y PkiPijVj = Vk + Y PkiPikVk + E E PkiPijVj ■ (40) 

3=1 i^kj i^k j+k i^k,j 

- . ' v » ' 

Ti: case ofj=k Ti\ case oij^k 
We point out that the term T± in the above equation recovers the desired signal lost in 
the 2nd stage, and the term T 2 removes the interference terms generated in the 2nd stage. 
Substituting (j5J) and fl5j) in PHI) , we can write 



(3) 

Vk 



K K 



K 



Xk 1 ~~ E PkjPjk - E Y pkjPjixi + n k ~Y^ Pkj 



n ■, 



j+k 



j¥=k l^j,k 



K 



+ E 

PkiPik 

i+k 
K 



K 



X k + Y Pkj X j + n k 
j+k 

K K 



3+k 

K K 

+ Y Y pkipij 

K 



K 



X 



+ Y pm + 



n i 



x k — Y, PkjPjk x k — E E PkjPjl x l + n k — Y, Pkj 



fl ; 



¥k 



K 



j¥=k j^k l=£j,k 

K K K 

+ E PkiPik x k + E PkiPik E P k j X j + E PkiPik n k 
ij^k i^=k jj^k ij^k 

K K K K K K K 

+ EE pkipij x j + Y Y p^pij Y pji x i + EE pmph 



n ; 



j+k i^k,j 



jjkk i^k,j 



¥j 



j+k i^k,j 



K 



K K 



K 



X k~Y PkjPjkXk ~YY PkjPjl x l + n k ~ Y Pkj 



j+k 



j¥=k ¥j,k 



j^k 







+ E 

PkiPik x k + E PkiPik Y Pkj X j + E PkiPi^k 
i^k iy^k j¥^k ¥k 



K K 



K K 



K 



K K 



K K 

+ E E p^Pij x j + Y Y PkjPjkXk + E E p^Pij Y pji x i + E E p^p^i ■ (4i) 

j+k i^k,j j^k i^k,j j^k i^k,j l^k,j j^ki^k,j 
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Appendix B: Expression for the mth stage filter output, m > 3 

From (jHJ), we can write y k m ^ for m > 3 as 

k ( K K K 

E Pkki _ E E Pkk 2 Pk 2 ki ~ E E Pkk 3 Pk 3 k 2 Pk2k 1 H 

ki^k ki=l \k 2 ^k,ki k 2 ^=ki k 3 ^k,k 2 

K K K K \ \ 

+ (- 1 ) m_1 H £ X! Y Pkkm-lPkm-lk^-i ■ ■ ■ Pk 3 k 2 Pk 2 k! ))Vk?- ( 42 ) 

&2^fcl fc37^2 fc m _2^fc m _3 km-l^k,km-2 J J 

Equivalently, (142]) can be written as 

2/i m) = yi m ~ 1} + (-i) m+1 E X X ••• E 

fel=l fe^fcl fe^fa fem-2^fcm-3 
' Pk 3 k 2 Pk 2 ki Vkx ■ 

(43) 

&m — 1 7^ ^ s & m — 2 



The first summation in the 2nd term on the RHS of (143]) can be split into two terms, one 
for k\ — k and another for k\ ^ k, as 

K K K K 

V™ = V™ + (-l) m+1 E E "' E E Pkk m -xPk m -ik m - 2 - • ■ Pk 3 k 2 Pk 2 kVk 

k 2 ^kk 3 =ik 2 k m ~ 2 =£k m -3 k m -ijtk,k m - 2 



T 3 

K K K K K 

+ (-i) m+1 E E E ■•• E E ■ ■ ■ Pk 3 k 2 pk 2 k x y kl ■ (44) 

klJ tk k 2 ^ki k 3 ^k 2 fc m _2^fc m _3 fc m _i^fe,fc m _2 



T 4 

For m > 4, 2/^ can be written in an alternate form as 

K K K K 

Vk = x k + (-i) m E E ■" E E Pkk m -iPk m -ik m - 2 • ■ ■ Pk 3 k 2 pk 2 kXk 

k 2 ^kk 3 ^k 2 km-^km-3 km-l^k,k m - 2 
" * ' 

T5: desired signal loss/gain at the (m-l)th stage output 
K K K K K 

+ (-!) m E E E ••• E E 

Pkk 7n —\Pk rn —\k rn — 2 Pk 3 k 2 Pk 2 k\ X k\ 
ki=ik k 2 ^k\ k 3 jtk 2 k m - 2 =£k m - 3 k m -i^k,k m - 2 
v ' 

T 6 : new interference terms generated at the (m-l)th stage output 
K K 

+ n k ~~ E Phki^ki + E E Pkk 2 Pk 2 k! n k! — • • • 
kiy^k ki=l k 2 ^k,k\ 

K K K K K 

+ (-1)™ E EE"' E E Pkkm- 2 Pk m - 2 km- 3 ' ■ " Pk 3 k 2 Pk 2 k x n kx . (45) 

fci=l fe^fci fc3^fc2 fc m _3^fc m _4 k m - 2 ^k,k m - 3 

Comparing the output terms in stages m—1 and m in Eqns. f|45l) and (JHJ), respectively, we 
can observe the following. 
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1. The desired signal loss/gain that occurred in the (m — l)th stage (i.e., T 5 in (|45l) ) is 
recovered/removed in the mth stage (see T 3 in (1441) and note that yjp in it has Xf.). 
In the process, new interference terms proportional to p m and additional noise terms 
proportional to p m ~ l get generated, due to all terms other than Xk in in T 3 . Note 
that for the case of m = 3, these interference and noise terms generated are given by 
| b, | and | b n | in Eqn. ( 14T1) in Appendix A. 

2. The new interference terms that were generated in the (m — l)th stage (i.e., Tq in 
(145]) ) are removed in the mth stage (see T 4 in (|44l) and note that y^ in it has x^J. 
In the process, new interference terms proportional to p m and additional noise terms 
proportional to p m ~ l get generated, due to all terms other than x^ in y^ in T 4 . Note 
that for the case of m = 3, these interference and noise terms are given by | E t | and | e n | 
in Eqn. (HTl) in Appendix A. 
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Figure 1: Average output SINR as a function of weight for the proposed G pw ^ m) filter. 
M = 1 (i.e., single carrier DS-CDMA), K = 20, P = 64, average SNR = 20 dB. No near-far 
condition. 




Figure 2: BER performance of various linear matrix filters - i) conventional filter Q,^ m \ and 
ii) proposed filter G p (m) . M = 1, K = 20, P = 64, average SNR = 15 dB. Near-far as well 
as no near-far conditions. 
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K = 20, P = 64, average SNR = 15 dB 
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Figure 3: BER performance of various linear matrix filters - i) MMSE converging filter 
G M (m) , and ii) modified MMSE converging filter, G PM (m) . M = 1, K = 20, P = 64, average 
SNR = 15 dB. Near-far as well as no near-far conditions. 



K = 20, P= 64 
average SNR = 1 5 dB 
No near-far effect 




Stage index, m 

Figure 4: BER performance of various linear matrix filters - i) filter, ii) proposed G p ^ 
filter, Hi) MMSE converging G^ (m) filter, iv) modified MMSE converging filter, G p/ /"% v) 
proposed weighted filter, G pw ^, and vi) conventional weighted LPIC filter in [9], G w ' m '. 
M — 1, K = 20, P — 64, average SNR =15 dB. No near-far condition. 
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Figure 5: Type-I Receiver for multicarrier DS-CDMA, M > 2. LPIC is done first on each 
subcarrier. Multicarrier Combining is done next. 
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Figure 6: Type-II Receiver for MC DS-CDMA, M > 2. LPIC is done after multicarrier 
combining. 
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Figure 7: BER performance of the proposed Type-I and Type-II receivers as a function of 
stage index m for MC DS-CDMA with - i) Type-I with G< m > filter, ii) Type-I with G p (m) 
filter, Hi) Type-II with G c ^ filter, and iv) Type-II with G p c(m) filter. M = 4, K = 20, 
P = 64, average SNR = 14 dB. Near-far effect as described in Sec. 13.61 
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